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Abstract 

The renormalization group equations (RGEs) of the mass matrices of quarks 

and leptons in a SO (10) model with two Higgs scalars in the Yukawa coupling 

are studied. This model is the minimal model of SUSY and non-SUSY SO(10) 

GUT which can reproduce all the experimental data. Non-SUSY SO(10) 

GUT model has the intermediate energy phase, Pati-Salam phase, and passes 

through the symmetry breaking pattern, 50(10) — > SU(2)l x SU(2)r x 

SU(4) C -» SU(2) L x U(1) Y x SU(3) C . Though minimal, it has, after the 

Pati-Salam phase, four Higgs doublets in Yukawa interactions. We consider 

the RGE's of the Yukawa coupling constants of quarks and charged leptons 

and of the coupling constants of the dimension five operators of neutrinos 

corresponding to the above symmetry breaking pattern. The scalar quartic 

interactions are also incorporated. 
PACS number(s): 12.15.Ff, 12.10.-g, 12.60.-i 
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Typeset using REVTfjX 



In these decades many informations on the quark lepton mass matrices have been accu- 
mulated. We are confronted with the very era when we should seriously consider a realistic 
model in the scheme of grand unified theories. On this ocasion, neutrino masses may be 
the window to the grand unified theories via heavy right-handed neutrinos. Along this sight 
we considered in [|T| and [Q] the SO(10) model where two Higgs scalars participate in the 
Yukawa coupling. Our SO(10) model with two Higgs scalar, 10 and 126, may be the first re- 
alistic model which successfully fit with the quark lepton mass spectra, Cabbibo-Kobayashi- 
Maskawa (CKM) ||, and Maki-Nakagawa-Sakata-Pontecorvo (MNSP) [§J mixing matrices 
based on the SO(10) framework. In these studies SO(10) invariant model is, of course, 
valid at the GUT scale and the data in our hand are those at much lower energy scale, at 
electro weak scale. So we must transport either of them to the other energy scale using the 
renormalization group equations (RGEs). In we bottomed up the electroweak scale data 
about six quark masses, three mixing angles and one CP-phase in the CKM matrix and 
three charged lepton masses to the GUT scale and, using the SO(10) GUT mass relations 
(see Eqs.(0) and flTSP) , obtained the neutrino mass matrix via seesaw mechanics J5| at GUT 
or at the intermediate scale. Then we toped down this mass matrix to the electroweak scale 
and checked whether it is consistent with the neutrino oscillation data and neutrinoless dou- 
ble beta decay. The results were very satisfactory in principle. However the theory is, of 
course, not conclusive. We adopted there some assumtions in order to extract the essential 
characters. One of them was the doublet-doblet splitting: that is, one of doublets is heavy 
relative to the other and our model in was the same in essence as the MSSM except for 
the additional mixing angle among heavy and light Higgs doublets. However, this doublet- 
doublet splitting is not sure to occur. For instance, in the same SO(10) GUT with two 
Higgs scalar, these two Higgs doublet may work in parallel to remedy the over- abundance 
of the leptogenesis 0. So in this paper we discarded this assumption and also incorporate 
scalar self coupling. In pay of this generalization we abandon in this paper the detailed 
data fitting and restrict ourselves in showing RGE's of quark-leptons mass matrices on more 
general ground. The renormalization of the neutrino mass operator was discussed by Babu- 
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Leung- Pant aleone and Chankowski-Pluciennik || for some simple models. Unfortunately it 
includes some calculational errors |J , though it does not affect our work, because we assume 
nuetrino couples to the Higgs doublets in a similar way to the up quarks. Our model is the 
most simple realistic model in SUSY and non-SUSY SO (10) GUT. At the GUT scale we 
have the Yukawa interactions to be given by 

£y = E (y^^ T B^H^ + ^Y^^B^i*^^^^^ + h.c^j , (1) 

where ty 1 (i = 1,2,3) is the 16-dimentional matter multiplet of the i-th generation, H and 
A are the Higgs multiplet of 10 and 126 representations under SO(10), respectively. B 
denotes the charge conjugation for SO (10) spinors : B = 7173757779. 
The Higgs scalar's quartic interactions with only H and A are given by 

V (if, A) = Ai {H»HJ (H»H U ) + iA 2 {H'HJ (a ww % W3W ) 
+ 1a 3 {H^ m5 ) (H V1 A—- **) 

(5!) 2 v ^/"iM2M3M4Msy y-^ L *-mv2vav4,vs ) ' a - u K^J 

Hereafter we consider non-SUSY SO(10) GUT explicitly. Gauge coupling unification needs 
the intermediate energy scale, Aj UlPf . Between the grand unification scale and the inter- 
mediate scale, the effective Yukawa interactions are given by 

- C Y = £ (YWFfQFi + Yg*>F?XFl + Y^F^F R + h.c) , (3) 

where F L and F R denote (2,1,4) and (1,2,4) in under G 22 4 = SU{2) L x SU{2) R x 
SU(4)a, respectively. And also, $, S and A R correspond to (2,2,1) in H, (2,2,15) and 
(1,3, 10) in A, respectively. Here we have assumed that suitably chosen U(1)h charge for- 
bids the Yukawa interactions like, Fl T $>F R and Fl T Y,F R , where 



(4) 



with 
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( A 

1 

And also the Higgs scalar's quartic interactions are given by 



V 



(5) 



V ($, E) = Ai [tr 2 + 2A 2 tr tr 



+ A, 



tr ( $S f 



tr 



4Xa 



tr fSS f 



(6) 



When the intermediate symmetry breaking occures, Ar have the vacuum expectation value 
and 



M. 



Ft 



(7) 



Below the intermediate scale, it includes four Higgs doublets : two doublets (0i and <p 3 ) 
come from (2,2,1) and (0 2 and 04 ) from (2,2,15). {0i, 2 } couple with up-type quarks 
and leptons, and {0 3 , 4 } do with down-type quarks and leptons |nj. In these multiple 
Higgs models we always encounter with the flavor changing neutral current (FCNC) problem 



as Glashaw and Weinberg remarked [ I2[ . However, by choosing the specific mass matrix 
textures, FCNC can be suppressed even in the case of multiple Higgs case jl3fl . Thus the 
low energy effective Yukawa interactions and the Higgs scalar's quartic interactions have the 
following forms, 



I.J 



(126) ' 



+ Y^f ¥ L <f> 3 e R + Y^f T L fae R + h.c. 



'■■J 



vtt)=\i E Kta (<f>%) (tin 

4 - a,6.c.d=l,2,3,4 

Here qi and £l are the left-handed quark and lepton doublets 



(8) 



(9) 



\' h J 



Also, 4> a are the Higgs doublets : 



( AO* \ 



From the hermiticity of @, the self coupling constants X a bcd satisfy 

^abcd = ^cdab = ^badc- 



(10) 



(11) 



(12) 



Eq.(^|) forbids the Yukawa and the scalar quartic interactions in Eqs.(|8]) and @ which are 
not invariant under the Z2 symmetry, 



^3, 94 



(13) 



and corresponding transformations in the fermion sector. When the electroweak symmetry 
breaking occures, (p a (a = 1,2,3,4) have the vacuum expectation values(VEVs), 



V2 



(a = 1,2,3,4) 



and 



M - -^y( 10 ) 4- _^2_y(i26) 
" ~ V2 u V2 u ' 

M _ V * v (10) V 4 (126) 



v 3 



y(10) + _^_y(126) 



V2 e V2' 

The Yukawa couplings above the intermediate energy scale are unified to 

1 y(io) = y(10) 

1 y(126) = ly(126) _ y(126) 

4v/2 F " 4 K 



(14) 



(15) 



(16) 



at the GUT scale. Here the above numerical factors in Eq.flTBD are necessarily to translate 
the SO (10) language into the G 2 24 language, see for example |TJ]]. And also, at the the 
intermediate scale, the Yukawa couplings above the electroweak scale are unified to 

y(10) = y(10) = y(10) = y(10) = yj,10) 

y(126) _ y(126) _ _ly(126) = _Iy(126) = yj.126) /-^ 

3 3 



and Eq.(^) is reduced to [[15 



M u = c 10 M<W + c 126 M( 126 ), M d = + M< 126 \ 

M D = c 10 M^ -3c 126 M (12e \ M e = M( 10 )-3M( 126 ). (18) 



Here 



Cio = f*/f 3 , C126 = , (19) 



and 



M (io) = ^yao) ; M ci*) = j* y a»>. (20) 

It goes from Eq.(^) that the Higgs quartic coupling constants are unified to 

^abcd 



4! 

^abcd 
Aaftcci 
A a ftcd 



Ai for a,b,c,d= 1,3 only. 
2A 2 for a, 6 = 1,3 and c,d = 2,4. (21) 
A 3 for {a, b} and {c, d} take any sets of {1, 2}, {1,4}, {3,2}, {3,4}. 
4A4 for a, b, c, d = 2, 4 only 



at the intermediate scale. Here the coupling constants which does not satisfy Eqs.(0) are 
vanished and excluded from the above rule. It should be emphasized that this model can 
be compatible with the large angle of atmospheric neutrino oscillation as far as we do not 
adopt any simplification not allowed in SO(10) framework [0]. For Ai 3> /i, Mr decouples 
and we must treat the four point interaction, 
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£ w = \ E E4 fl,6) N(^3) + /l ' c - 

a,6=l,2 ij 
a, 6=1, 2 ij 

| + 3>{?) + } + h.c, (22) 



where i,j are flavour indicies. The reason why a, b run over 1,2 is that it comes from the 
Dirac neutrinos (the third line in the right-hand side of Eq.(5)). Thus the effective light 
neutrino mass matrix is given by 

M„ = ~ E ^%lvl (23) 

a,&=l,2 

We have two Higgs doublets in each up quarks, down quarks, Dirac neutrinos and charged 
leptons and k^'^ (a, b = 1,2) in ours are the generaslization of in m. We have no 
term corresponding to k( u ) in [HI. 

In the non SUSY SO(10) with two Higgs scalars, the one loop RGEs for the effective Yukawa 
couplings first at the energy region between the grand unification scale and the intermediate 
sale are given by : 

igtt 2 ^^ = (V( 1( M 10)t + -y£ 26) y^) y( 10) 



i 



, y(10) fy(10)y(10)t , ^ /y(126)y(126)f , y (126) y (126)f\ 

'i^j^i* 1 F ' ^ y F F ' R R J 

+ 4 tr (rrrr 1 ) rP + (% + \ot R + f afc) if ', (24) 



2 rfy F _ f v (10) v (10)t i 15 V (126) v (126)t\ v (126) 

y(126) fy(10)y(10)t , 1^ /y(126)y(126)f y(126)y (126)f\ 1 
~ .F 1 I F F 1 ' ^ y F F ' R R, J ( 

+ te (r^M 1 26)t ) if 26 ' + (!& + 1«5» + ^Ac) rr>, (25) 

J V (126) , ,r 

)"•* R I -. r(1(l), rfmlt . 10 /, ,(126l,,fl26lt . T ,d26W,d2fiH\ 1 -^(126) 



16n2 dF~ = { Y ? 0)Y ? m + f ( y i 126)y i 126)t + 4 126) 4 126)t ) } ^ 



y(126) I y(10)y(10)| , 15 /y(126)y(126)f y (126)y (126)f\ 

i? I i* 1 F ' a \ F F ' R R J 



4 

+ tr (Y^'Y^) Yr + (§& + if, (26) 
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where g 2 L, 92R and g^c are the SU(2) L , SU(2) R and SU(A) C gauge coupling constants, 
respectively. At the second stage, the energy region between the intermediate sale and the 
weak scale, the one loop RGEs for the effective Yukawa couplings are given by : 

16n 2 ^— = 3tr(Yj 10 )yj 10 )t)y u ( 10 ) + 3tr(Fj 10 )Yj 126 )t)Y u ( 126 ) 



dt 



+ 2 



I (y(10)y(10)f + y(126)y(126)f 

+ Y] 10) Yj 10)t + yj 126 )yj 126 )t^ y(io) 

+ ^ (10) (Vj 10)t n (10) + y^V ] ) , (27) 

i67r 2^_ = 3 tr(Fj 126 Vj 126 )t)Yj 126 ) +3tr(Fi 126 )Yi 10 )t)F u ( 10 ) 



dt 



1 /V(10)y(10)t , y(126)y(126)f 

2 ^ u u u u 

+ F] 10) Fj 10)t + yj 126 )yj 126 )t^ y(i26) 

+ Fi 126 ) (Yi 10)t Yi 10) + Y^Y^) , (28) 



+ 2 



,(10) 



1&r dT = { 3tr(y « ° hfm] + tr ( y '. (10,y '. (10,t )} y i 10 ' 
+ {3tv(Yfhf m ) + tr(y;<'"iy;<' 2 «>')} rj 1261 



+ 2 



1 (V(10)y(10)t , y(126)y(126)t 

2 ^ « it ~ u u 

, v (10) v (10)t , v (126) v (126)f\ v (10) 

+ ^ 1 °H^ 1 °M i nri 126 %r ) ), (29) 



^v( 126 ) 

167r 2^_ = / 3 tr(Yi 126 Vi 126)t ) +tr(Y e (126) Y e (126)t )} Yi 126 



+ {3tr(y;< 126) yj 10)t ) + tr(y e ( 126 )y e ('°)t)} yj 10) 

+ I (y(10)y(10)t + y(126)y(126)f 



2 

, v (10) v (10)f , v (126) v (126)f\ v (126) 



+ Y W( Yi m Yi io )+Y w Y OX>) i (go) 

lQ^^jp- = {3tr(ri 10 Vi 10) Vtr(r e (10 Vi 10 )t)}r e ( 10 ) 
+ {3tr(ri 10 Vi 126) Vtr(ri 10) ri 126)t )}r e (126) 

+ I fy(10)y(10)t + y(126)y(126)f\ y(10) 
e e e y e 

+ rj 10 > (y e (i°)ty e (io) + y e (i26) tye (i26)^ ) (31) 

1 67T 2^£!!! = { 3tr( y^6) r jl26) t) +tr (y e (126) n (126) t) j ye (126) 

+ {3tr(ri 126 Vj 10)t ) + tr(ri 126 )ri 10 )t)}y e ( 10 ) 

- + ™&) n {126) 

+ I (y(10)y(10)f + y(126)y(126)f\ y (126) 
2 V e e e e / e 

+ Y" e (126) (y e ( 10 )ty e ( 10 ) + y e (i26)tyji26)^j _ (32) 
In Eqs. (f27l) -(32l), the first two terms (before gf terms) are the contribution of fermion loop, 



the third term that of gauge loop, and the remaining two terms those of Higgs loop. The 
second term indicates the mixing of two Higgs doublets in the fermion loop correction. 
These formulas, of course, are reduced to one-loop RGEs for the standard model if we set 
either y( 10 ) or y( 126 ) zero fl6fl . The one loop RGEs for the scalar quartic couplings in our 



model are [17 



167T ^ ] (2A a fe mn Xnmcd Xabmn^cmnd ^amnb^mncd 

^amnd^cnmb Xamcn^mbrid) 3 ^3(^2 9y^) ^abcd 

+ 9 (3^2 + 9y) $ab8cd + 36 g\ g Y ^S ad 5 bc - ^S ab 5 cd 

^ ] i^mbcd-A-am ^amcd-^mb ^abrnd-^-cm Xabcm-^-md) 
m,n=l, 2,3,4 

-48H abcd {a,b,c,d= 1,2,3,4), (33) 



Here 



A ab = tr (3Y?Y? + 3Y a d % d + Y?Y*) , (34) 



and 



H abcd = tr (3Y d ^Y c u Y b u X u + 3Y d % d Y d % d + Y«%*Y?Y* 

+ ZY?Y?Y«Y* + 3Y d % d Y c u % u - 3YfY d Y b U % u ~ 3Y a u % u Y b ^Y d ) , (35) 



with 

V u _ y-(10) V u _ v (126) v d _ y(lO) yd _ y(126) 

Kf = y e ( 10 ), y 4 e = y e ( 126 ), and otherwise zero. (36) 

As for the RGEs of VEV's of Higgs fields, there may be some conflicts. Someone consider 
them constants || [|T^] [I5|, and other ones make them evolve as dy ^ ^ for the simplest 



case |[20|| . The situation may depend on what scale and what object we consider. Here we 
adopt the standpoint that the RGEs of v a are those of 0° [|16j |22|. That is, in our case, 

16vr 2 ^ = -3tr (Y^Y^) - 3tr (Y^Y^) v 2 



IQ^^l = _ 3tr (y u (126) y a26) t ) ^ _ 3tr (yjl^yjlOt 



+ ('t02 + 7^K> ( 37 ) 



4 . 



tit 



«1 



+ ( 7#I + -*9 2 y ) v 2 , (38) 



16vr 2 ^ = -tr (3ri 10) Vi 10) + r^ty^o)) U3 



tr (3Fi 10) Vi 126) + Y^Y^ 
9 2 3 

-09 + - 



+ (t^ + 7^) v 3, (39) 



167r 2 ^i = -tr f 3 yj 126 )tyj 126 ) + y e (i26)t Fe (i26)\ ^ 



- tr (3Fi 126)t Fi 10) + yO»)ty e (iO)) ^ 



+ (~<? 2 2 + j<&) «*■ (40) 
Here we have replaced the RGEs of <fi a with v a . It goes from Eq.flTo]) and Eqs.(P7|)-(PDD that 



167r2 ^i = _JL tr (yao)y(i 2 6)t _ y^yOOM) (, 2 yao) 



y(i26) 
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_i_ - (yj 10) y u (10)t + yj 126) y u (126)t 



+ Y } w) Y} m + Y} 126) Yj 12 ^)M u 

+ M u (Fi 10)t Fi 10) + Y^Y^) , (41) 



2 dM d / 2 1 2 \ 
16vr — = - [8g 3 - g^J M d 



+ I M 10 )yJ 10 )^ + yjl26)yjl26)f 

+ ^ 10 Vi 10)t + ri 126 Vi 126)t )M, 

+ M d (yj 10 ^ + yj 126)t y d 126 ) , (42) 



9 dM e 9 , , 

16tt 2 — - = -3glM e 
dt yy 



+ i (y ao, y ao„ +yr)y « t)Me 



+ M e (yP°)tyW + y;<r26)tyUM>) . (43) 

It should be remarked that RGEs destroy the transpose-invariance of mass matrix possesed 
at GUT. At the region Aj ^> we must treat the RGEs of k of Eq. (|22|) , giving the following 
forms; 

16tt 2 ^ = etr (ywywt) 

+ 3tr (yi 126) y M (10)t ) (/t (L2) + « (2 - 1} ) 

+ i (A mi «W) + A m2 K (1 ' 2) + A 1212 ^ 2 < 2 >) 

+ k^' 1 ) (y e ( 10 )y e ( 10 )t + Fe (i26) Fe (i26)t) T j ; (44) 

,7^(2,2) 

i67r 2^ = Q tr / r (126) r (126) t A (2,2) 

dt \ u u ) 



+ 3tr (yj 10 )yj 126 )t) (V 1 - 2 ) + k( 2 ^ 

- 3^ 2 V 2 ' 2) 

+ I (A 2222 k;( 2 ' 2 ) + A 2 i2i/« (1 ' 1) ) 

+ 1 |(Y" e ( 10 )Y" e ( 10 )t _|_ ya26)yjl26)t^ R (2,2) 
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+ (ri 10 Vi 10 )t + y e ( 126 )y e ( 126 )t) r j , (45) 
i 6vr 2^!! = 3tr (yao )y ao) t + Y^Yi 12 ^) 

(JjL 

+ 3tr (rj 10) yj 126)t ) «: (L1) 
+ 3tr (rj 126) yj 10)t ) «: (2 - 2) 
- y 2 (2^^) + K w) 

+ i (A 1122 «: (1 ' 2) + A 1221 ^') 

+ I |(y e ( 10 )y e ( 10 )t + ya26)y e (126)t^ ^(1,2) 

+ ^' 2 ) (rj 10 )yi 10 )t + r e ( 126 )y e ( 126 )t) T | , (46) 

levr 2 ^— = 3tr (ri 10) ri 10 )t + yj 126 Vj 126 )t) 
+ 3tr (rj 10 )yj 126 )t) 
+ 3tr (rj 126) yj 10)t ) «: (2 - 2) 

-^(^ 2 ) + 2^ 2 ' 1 )) 

+ i + A 2112 ^ (1 ' 2) ) 

+ I |^y e (10)yj:iO)t + yj[126)yj[126)t^ ^(2,1) 

+ k^ 1 ) (Fi 10 Vi 10 )t + r e ( 126 )y e ( 126 )t) T | , (47) 

The different factor 2 of the coefficients of g\ in «( 1,2 ) and k^ 2 ' 1 ^ comes from the different 
contribution of W-boson and Z-boson loop correction on k^ 1 ' 2 ) and k^ 2 ^. The self coupling 
contributions come from the diagram Fig.l. 

Substituting the above equations into Eq.(18), we obtain the RGE of the light neutrino mass 
matrix: 

16vr 2 ^ = i{3(, 2 + , 2 )M, 

a,b,c,d=l 1 2 

+ (rj 10 Vj 10 )t + r e (i26)y e (i 2 6)t) Mu 

+ m v (y^ywt + ri 126 )ri 126 )t) T | . (48) 
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Thus we have obtained RGEs of SO(10) GUT with two Higgs scalars. These formulas are 
especially important for essentially multiple Higgs models where two Higgs doblets in up or 
down type quark-lepton mass matrices can not suffer doublet-doublet splittings. Even in a 
model where the doublet-doublet splitting occurs, our results are useful for estimating the 
threhold effects due to the heavy pairs of Higgs doublets. We calculated RGE's explicitly 
in non-SUSY SO(10). Those of SUSY SO(10) are obtained analogously by considering the 
contributions of the superpartners of non-SUSY contents except for the scalar self couplings. 
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FIGURES 



VL <t>c 




FIG. 1. The Higgs scalar self couplings contribution to the RGE's of k. 
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